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I. INTRODUCTION 

Quantum interference phenomena are essential when 
developing mesoscale electronic devices. Quantum con- 
fined geometries conceived for such studies may consist 
of two-path interferometers, 1,2 parallel quantum dots, 3 
coupled quantum wires, 4 side-coupled quantum dots, 5,6 
or Rashba double dots in a ring. 7 These coupled meso- 
scopic systems have captured recent interest due to their 
potential applications in electronic spectroscopy tools 8 
and quantum information processing. 9 Nevertheless, a 
study of microscopic magneotransport behavior of the 
transient current flow in an interacting window-coupled 
double quantum wire system is still lacking. 

In the presence of a magnetic field perpendicular to the 
plane of the wires, the energy spectra have been stud- 
ied pointing out the complex structure of the evanescent 
states of the system in homogeneous 10 and inhomoge- 
neous 11 double wires (DW). It was shown that the step- 
wise conductance increasing and decreasing features can 
be changed by the applied magnetic field and the height 
of the barrier between the wires. 12 Moreover, the dynam- 
ics of the transfer processes for single-energy electron 
spectroscopy in coupled quantum states has been con- 
sidered with window coupling potential experimentally 13 
and theoretically. 14 

In a closed time-dependently driven quantum system, 
the Jarzynski relation may be derived without quan- 
tum corrections by introducing the free-energy difference 
of the system between the initial and final equilibrium 
state. 15,16 When the system is coupled to the reservoirs, 
the Jarzynski relation can be derived using a master 
equation approach. 17,18 Different approaches were pro- 
posed based on the quantum master equation (QME) to 
study interaction transport effects. 19 21 The time evolu- 
tion of the system described by the QME consists of two 
parts: The Hamiltonian describing the system induces a 
unitary evolution of the reduced density matrix, and the 
dissipative part describing the properties of the environ- 
ment or reservoirs. 22 

To study the time-dependent transport properties, the 



assumption of Markovian dynamics and rotating wave 
approximation lead to different types of master equa- 
tions of the density matrix for the study of steady-state 
currents by neglecting memory effects in the system, 23 
in which the diagonal and off-diagonal elements of the 
reduced density operator are decoupled 24 or assuming 
an infinite bias regime. 25 However, the transient time- 
dependent transport, which carries the coherence and 
relaxation dynamics, cannot be generally described in 
the Markovian limit. An accurate numerical method 
for the nonequilibrium time-dependent transport in the 
interacting nanostructures is desirable, which can ver- 
ify various approximation approaches. A non-Mar kovian 
density-matrix formalism involving the coupled elements 
should be considered based on the generalized QME 
(GQME). 26 30 It has been confirmed that the Marko- 
vian limit not only neglects the coherent oscillations, but 
also the rate at which the steady state under this limit 
significantly differs from the non-Mar kovian results. 30 

In this work, we investigate how the interplay of the 
magnetic field and the electron-electron (e-e) interaction 
affects the quantum interference of the parallel quantum 
wires through a coupling window with a time-dependent 
switching-on coupling to the leads. The central finite 
DW system is connected to semi-infinite leads of the 
same width. To explore the switching-on time-dependent 
transport behavior through the sandwiched DW system, 
we shall explicitly construct a transfer Hamiltonian that 
is spatially located at the system-lead contacts and with 
a certain distribution in the energy domain. Due to the 
finite size of the DW system, the Coulomb correlation 
could play important role in the transport. Appropri- 
ately tuning the above physical parameters, we obtain the 
transient as well as the quasi-steady state electric current 
using a non-Markovian GQME method. This allows us to 
explore quantum interference features of the dynamical 
transient currents through the tunable window-coupled 
DW system. 

The paper is organized as follows: In Sec. II, we present 
the model describing the window-coupled DW system 
based on the GQME theory. Section III presents our 
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numerical results and physical discussion, 
remarks are addressed in Sec. IV. 



Concluding the length of the DW system and V c (y) = |m*^Q?/ 2 is a 
parabolic confining potential. It is convenient to rewrite 
the non-perturbed single-electron central system Hamil- 
tonian as 



II. MODEL AND THEORY 

Quantum transport in an open system acted upon 
by a time-dependent potential has been consid- 
ered in different systems such as time-dependent 
quasibound-state features, 31 ' 32 quantum pump in 
Luttinger liquids, 33 photon-associated transport in 
nanostructures, 34 36 the Kondo effect in a double 
quantum dot— quantum wire coupled system, 37 ac-field 
control of spin current, 38,39 and transient current dy- 
namics in nanoscale junctions. 40,41 The rapid progress 
of nanoelectronics and information technologies has 
prompted intense interest in exploiting the quantum 
interference transport properties of correlated electrons, 
in which the coupling between the mesoscopic subsystem 
could be manipulated by an applied external magnetic 
field. Furthermore, the increasing interest in fast dy- 
namics in mesoscale systems and time-resolved detection 
of electrons via a nearby detector strongly motivates 
investigations of interacting time-dependent transport. 
It is thus warranted to explore the magnetotransport in 
a central system that is weakly-coupled to the leads by 
switching-on time-dependent potentials located at the 
system-lead junctions. 
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for defining the effective cyclotron frequency 
uj 2 in terms of the two-dimensional cyclotron frequency 
uj c = eB/m* . The typical length scales of the system 
along the x and y directions are characterized by the 
two-dimensional magnetic length I = (h/m^Uc) 1 ^ 2 and 
the modified magnetic length a w = (ft/m*^) 1 / 2 respec- 
tively. 

Utilizing the microscopic single-electron eigenfunctions 
of the system ^n( r ) allows us to express the system 
Hamiltonian in the spectral representation 42 



n 1 1 
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where E n stands for the eigenvalues of the central system 
and the dummy index n refers to the quantum numbers 



Considering the parabolically confined semi- 



infinite leads, one obtains the single-electron Hamiltonian 



Yl / d( i e n y (QM l n y , q )(^n y , 



(6) 



A. Single-electron Model 

One starts from an open quantum system described by 
a single-electron time-dependent Hamiltonian 



h(t) =h + h T (t). 
Therein, the first term 



h = hs^ riJ 

l=L,R 



(i) 



(2) 



indicates a disconnected single-electron Hamiltonian de- 
scribing the central system by h$ and the biased leads by 
hi with I referring to the left (L) and right (R) leads; and 
the second term h^(t) stands for a switching-on time- 
dependent transfer Hamiltonian connecting the central 
system and the leads. The h$ contains a disconnected 
Hamiltonian ho and an envelop potential Vow( r ) de- 
scribing the embedded double quantum wire subsystem, 
namely 



• Vbw(r). 



(3) 



Here h$ = p 2 /2m* + V con f(x, y) is composed of a kinetic 
term with canonical momentum p = p + eA with vec- 
tor potential A = (0, —By, 0) and a confining potential 
Konf(^2/) = V c (x) + V c (y), where V c (x) denotes a hard- 
wall confining potential at x = ±L x /2 with L x being 



in which q stands for the continuous wave number along 
the transport direction and n l y denotes the transverse 
subband index with I referring to either of the two leads. 
We assume the contact is gradually switched on in time 
and calculate the time-dependent reduced density oper- 
ator of the sample using the GQME. The DW system is 
coupled to the leads by introducing the off-diagonal time- 
dependent transfer Hamiltonian /it(^) = hrfi(t) + h^(t), 
where 



h l T (t) 



fdq X l (t) (T< 



h.c.) 



(7) 



with T l qn being the coefficients connecting the eigenstates 



in the system ip^ and the leads ip q . Explicitly, we express 
the switching-on contact function in the I lead as 



x l (t) = 6(t-t ) 



e t(t-t ) _|_ i 



(8) 



such that the coupling between the central DW system 
and the leads is switched on at t = to and the parame- 
ter 7 indicates the switching rate of the coupling. The 
current will flow through the system once the switching- 
on contacts between the device and the leads have been 
established. 



B. Many-electron Model 

The Coulomb interacting many-electron states (MES) 
of the isolated sample are derived with the exact diago- 
nalization method. 43 The chemical potentials of the two 
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leads create a bias window which determines which MES 
are relevant to the charging and discharging of the sam- 
ple and to the currents, during the transient or steady 
states. The many-electron Hamiltonian 



H(t) = H + H T (t) 



(9) 



consists of a disconnected many-electron system Hamil- 
tonian 



l=L,R 



(10) 



and a time-dependent transfer Hamiltonian H^(t). The 
central system Hailtonian H$ = H$ + H$ contains a 
kinetic term H$ = J2 n E n S n d n with discrete single- 
electron energies E n and a Coulomb interaction term 



n\m' ;n,md n / d m / d n d m , (H) 



n' ,m' n,m 



where we have introduced the electron creation (annihi- 
lation) operators in the system d\ (d n ). The two-electron 
matrix elements 

Vn' ^rn' ;n,m (1^) 

dvdv'^ n ,(vY^ m ,{v'yV{v - r')^(r)^(r'), 



expressed by the single-electron state (SES) basis, are 
derived for the Coulomb interaction potential 



V(r - r') 



^7re e r y/( x - x') 2 + (y - y') 2 + rj 1 



(13) 



with e r and r] being, respectively, the relative dielec- 
tric constant of the material and the infinitesimal con- 
vergence parameter. Below we define the dummy index 
q = {n l y1 q) and J dq = J2n y I for simplicity. The 
many-electron lead Hamiltonian can be expressed in the 
following form 

Hi=f dqe\q)c l q ] c l q . (14) 
The second term in Eq. (9) is expressed explicitly as 
H l T {t) = X \t)Y J [ dq [c l qT l qn d n + di{T l nq Yc l q 

(15) 

describing the transfer of electrons between SES of the 
the system \n) and the leads \q) through the coupling 
coefficients T^ n , given by 



L qn 



I drdr%(r'yg l qn (r,r'W s n (r). (16) 



Therein, the coupling function 

g l qn(ry) = g l exp [-S l x (x - x'f - S l y (y - y'f] 

x exp (-A l n (q)/A) (17) 



containing the system- lead SES energy spread A l n (q) = 
\E n — e l (q) | making the connection of any two SES at the 
contact region in the energy domain. 29 The spatial cou- 
pling range in the leads is governed by 5 l x and S l y . We have 
considered the energy interval [pr — A, /i^ + A] to define 
an active window in the energy domain A# = Ap + 2A 
that involves all the possible states in the central system 
that are relevant to the transport. It should be men- 
tioned that only the transverse part of the wave function 
in the semi-infinite leads is normalizable. To get rid of 
all length scales variation with magnetic field, one needs 
to fix g l aw 2 in units of energy and then calculate g l . 



C. GQME Formalism 

In this subsection, we formulate the time evolution of 
the MES when the system contains a number of elec- 
trons for the study of interacting time-dependent trans- 
port properties based on the GQME formalism. 44 To 
take into account the many electrons in the system, we 
construct a Fock space by selecting the number of the 
TVses lowest single-electron states and the Nmes — 2 NsES 
many-electron states within the active window A#. In 
the occupation representation basis, the noninteracting 
MES 



\a) 



L l 5 L 2 1 ' 



3 ) 



(18) 



contains the labels i% = 0, 1 indicating the occupation of 
the n-th SES of the isolated central system within the ac- 
tive window. The corresponding energy of the noninter- 
acting MES £ a = E n i% can be obtained by summing 
over the occupied SES. 

The time-evolution of the many electron system under 
investigation obeys the Liouville-von Neumann (quan- 
tum Liouville) equation 45 



dW(t) 
dt 



-j.[H(t),W(t)}, 



(19) 



where the full density operator W(t) can be operated 
upon by a projector to yield the reduced density oper- 
ator (RDO) by taking trace over the Fock space in the 
leads p(t) = Tr L Tr R W(t), with p(t ) = p S - 46 The initial 
condition W(t < to) = PlPrPs is in terms of the equi- 
librium RDO of the disconnected lead I with chemical 
potential /i/, given by 



-PiHt-fjHNt) 



pi 



Tri{e~^ H i-^ N i)} 



(20) 



with I referring to the L and the R leads. This allows 
us to find the equation of motion for the RDO of the 
following form 47 

^ = -i£ eS p(t) + jf dtX(t, t')p(t'), (21) 

where C e R stands for the effective Liouvillian and JC(t,t f ) 
denotes the integration kernel. 47 
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Using the exact diagonalization method, we diagonal- 
ize the interacting system Hamiltonian Hs in the MES 
basis of the noninteracting system in the Fock 

space. Since we are dealing with an open system with 
variable electron number, one has to include all sectors 
containing zero to Nses electrons. This yields a new in- 
teracting MES basis {\p)} with 



1^0 = ^U^Ot) 



(22) 



connected by the A/mes x A/mes unitary transformation 
matrix U^ a . A basis transformation of the interacting 
many-electron coupling matrix T l (q) = U^T l (q)U and 
the insertion of the diagonalized matrix representation 
of the interacting Hs allows us to obtain the RDO in 
the interacting MES basis p = W pU. Expressing the 
interacting many-electron coupling matrix T in the in- 
teracting MES 



r l (q) = J2 r U^)M 



(23) 



with T* u (q) = En T ng(/ i l4l^) in terms of the single- 
electron coupling matrix T^ q , one can obtain the trans- 
formed GQME 



dp(t) _ i 

— _--[If s ,p(f)] 



(24) 



-i E I d * ([r l (q)M l q(t)} +h.c). 

l=L,R J 

Here we have defined the effective interacting coupling 
operator 



£l l q(t) = Ul(t) f ds X l (s)U l q (s) 

Jtn 



x exp 



h 



(t-s)e\q) 



U s (t) (25) 



with 



n^(s) = u s (s) 



(V)V(s) [I"/' «q))} 



p(s)(T i y f(e(q)) Ul(s) 



in which Us(t) = e rlH s(t-t )/n d eri otes the time evolution 
operator and f l (e(q)) = {exp[e(g) — m] -f 1} _1 indicating 
the Fermi function in the I lead at t = to. In the numer- 
ical calculation we shall select to = for convenience. 

Taking the statistical average over the Fock space 
(Qs(t)) = Tr{W(t)Q s } of the charge operator Q s = 
e d) n d n in the coupled central system and using the 
identity p(t) = TrLTrR{W(t)}, one may express the sta- 
tistical averaged time-dependent charge as 



(26) 



This allows us to define the time-dependent net charge 
current flowing through the central DW system 

I Q {t) = d ^^=I L {t)-I R {t). (27) 

The charge current injected from the I lead to the system 
is given by 



n fi 



dt ' 



(28) 



in which we express the current in terms of the time 
derivative of the reduced density matrix elements in the 
interacting MES basis: 



d P l nu. 
dt 



x l (t) 



dq (p\[T l (q),n l q (t) 



h.c. 



It is straight forward to obtain the interacting many- 
electron charge distribution in the DW system 

Q(v,t) =e^C,(r)^(r)^p^(i)( M Mi,rf„|^). (29) 



Below we shall show our numerical results of the net time- 
dependent charge current Iq{t) through the central DW 
system. It is an algebraic sum of the left current (in- 
dicating the charge current from the left lead to the right 
lead) and the right current (indicating the charge 

current from the system to the right lead). We shed 
light on the transport dynamics by analyzing the time- 
dependent many-electron charge distribution Q(r,t) in 
real space. 



III. RESULTS AND DISCUSSION 

We numerically solved the GQME to investigate the 
dynamical time-dependent magnetotransport of elec- 
trons through a central finite system of length L x = 
300 nm with magnetic length I = (h/ \eB)) 1 ^ 2 = 
25.67[£?(T)] -1 / 2 nm. The central system is transversely 
confined by a parabolic potential with characteristic en- 
ergy HQ.Q = 1.0 meV. This supplies the modified magnetic 
length 



1/2 



ra*ft / V 1 + (e£/(ra*^ )) 2 
33.74 

nm, 



1/4 



{/1 + 2.982[£(T)] 2 



(30) 



and the typical width of the confined system for the low- 
est subband electron is L y « 67.5 nm. We assume GaAs 
parameters with electron effective mass m* = 0.067ra e 
and the background relative dielectric constant e r = 12.9. 

Figure 1 schematically illustrates the window-coupled 
DW system scaled by a w . The embedded DW system is 
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FIG. 1. (Color online) Schematic illustration of the potential 
defining the window-coupled DW system. hQo = 1.0 meV, 
B = 0T, a w = 33.74 nm. 

described by Vow( r ) = Vmb(2/)+Vcw(^ 5 v) that contains 
a middle barrier 

V M B(y) = V eM-Poy 2 ) (31) 

with Vq = 18.0 meV and /3q = 0.3 nm -1 , as well as a 
coupling window potential 

Vcw(x, y) = -V eM-PW - P 2 y y 2 )- (32) 

The coupling constant g l aw 2 = 60 meV, and the contact 
size parameter S l x = 5 l y = 4.4 x 10 -4 nm -2 

In the following calculations, the temperature of the 
reservoirs is fixed at T = 0.5 K, and the states within 
the bias window before switching-on the coupling are as- 
sumed to be unoccupied. The coupling between the DW 
system and the leads is characterized by the switching 
rate 7 = 1.0 ps -1 , and the nonlocal coupling strength 
is fixed as T l = 4g l aw 2 /(d^Sy) 1 ^ 2 = 54.5 meV-nm 2 . 
The bias voltage is fixed leading to a bias window 
eVbias = A/i = 0.9 meV, and the extension parameter 
A = 0.3 meV is selected referring to a window of rele- 
vant states Ae = A/i + 2 A = 1.5 meV. 

The energy spectrum of the leads as a function of wave 
number q scaled by a" 1 is shown in the left panel of 
Fig. 2. The bias window A/i is located in the first sub- 
band, whereas the extended active bias window covers 
the evanescent modes below the first subband and the 
threshold of the second subband. The energy spectrum 
of the window-coupled DW system as a function of the 
single electron number n is shown in the right panel of 
Fig. 2 containing five SESs in the window of relevant 
states A#; the three lowest states are in the bias win- 
dow A/i whereas the two highest states are in the upper 
extended window [/i/, , /i/, + A] . 

In Fig. 3, we show the time-dependent charge cur- 
rent for the case of magnetic field B = 0.5 T with and 
without e-e interaction, denoted by Iqj and Iq^ re- 
spectively. The noninteracting left and the right cur- 
rents are also presented for comparison, denoted by Il,q 
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FIG. 2. (Color online) Energy spectrum of the leads (solid 
red) versus wave number q (left panel); and energy spectrum 
of the window-coupled DW system (cross dot) versus the SES 
number n (right panel). Magnetic field B — 0.5 T, and the 
chemical potentials in the the leads are \il — 1.65 meV and 
Hr = 0.75 meV (dashed green) such that A/i = 0.9 meV. The 
window of relevant states Ae = 1.5 meV is defined by the 
dotted blue lines. 

and Jr 5 o respectively. We have selected j3 x = 0.02 
nm -1 and j3 y = f3o such that the length of the cou- 
pling window L w is 100 nm. In addition, the coupling 
constant is g l aw 2 = 60 meV and the contact size pa- 
rameters are 5 1. — Si — 4.4 x 10 -4 nm -2 such that 

x y 

the coupling strength T l = 54.5 meV-nm 2 and the ef- 
fective lengths of the system-lead coupling potential are 
L l c x = L l c y w 95 nm. Below we shall show that the 
time-dependent charge current manifests different trans- 
port mechanisms in the short-time and the long-time re- 
sponse. 

In the short-time response regime, shown in Fig. 3(a), 
the time-dependent charge current is increased and man- 
ifests rapid oscillation with period r s w 1.9 ps exhibiting 
quantum interference dominant features. In this regime, 
the noninteracting approach could be a good approxima- 
tion for analyzing the transient time-dependent transport 
properties. In this short-time regime, the interacting and 
the noninteracting currents are almost the same before 
20 ps with negligible right charge current implying effec- 
tive charging and quantum interferance dominant trans- 
port feature. The right charge current is significantly 
increased after 20 ps. At around t = 35 ps, the difference 
between the interacting and noninteracting currents be- 
comes 0.1 nA (the Coulomb correction is ~ 10%), and 
the right charge current is increased to 0.18 nA. 

In the long-time response shown in Fig. 3(b), the 
charge current displays slow quasi-periodic oscillation 
with period 77 w 39 ps approaching a steady current. 
The slow oscillation behavior in the time-dependent cur- 
rent implies that the quantum interference feature is sup- 
pressed whereas the Coulomb interaction effect is en- 
hanced. At time t = 250 ps, the interacting steady cur- 
rent (~ 0.17 nA) is much higher than the noninteracting 
steady current (~ 0.015 nA). The mean charge of the DW 
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FIG. 3. (Color online) The interacting net current Iqj (solid 
red), the noninteracting net current 7q,o (solid blue), the non- 
interacting left current Il,o (dashed green), and the noninter- 
acting right current 7r,o (dashed light-blue) are plotted as 
a function of time: (a) short-time response; (b) long-time 
response. B = 0.5 T, L w = 100 nm, the coupling con- 
stant g l aw 2 = 60 meV, and the contact size parameters 
5 l x = S l y = 4.4 x 10" 4 nm~ 2 . 



system is monotonically increased in time (not shown), 48 
and the mean charge of the interacting MES (~ 0.8e) is 
approximately twice that of the steady mean charge of 
the noninteracting MES (~ 0.4e). This indicates that the 
empty-state initial condition ensures that the Coulomb 
interaction facilitates to drag the electron dwelling in the 
DW system through the window of relevant states, and 
thus enhances the steady current. 

In order to get better understanding on the transient 
dynamical transport, we present the spatial distribution 
of the many-electron charge at t = 10, 25.6, 100, and 
200 ps in Fig. 4, labeled by a-d in Fig. 3, respectively. 
When the system-lead coupling is switched-on with for- 
ward bias, the electrons are incident from the left lead 
into the system with transversely symmetric distribution 
(not shown). At around t c± 10 ps, the electrons located 
in the lower wire favorite to make inter-wire backward 
scattering to the upper wire, exhibiting a fully quan- 
tum mechanical feature. Later on, the electrons perform 
an opposite inter- wire backward scattering feature to the 
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FIG. 4. (Color online) The many-electron charge density for 
the noninteracting (left panel) and interacting system (right 
panel) for B = 0.5 T. The other parameters are the same as 
Fig. 3. 



lower wire at t ~ 25 ps, and this feature is only slightly 
enhanced by the Coulomb interaction. However, in the 
long-time response regime, the inter-wire scattering for- 
ward and backward effects are both enhanced. At around 
t = 100 — 200 ps, the noninteracting window-coupled DW 
forms a quasi-isolated four cavities, the window coupling 
effect is significantly enhanced by the Coulomb interac- 
tion. It is interesting that the electron can form a quasi- 
bound state in the coupling window at t ~ 200 ps. When 
the DW system approaches steady-state transport in the 
long-time response regime, the total charge in the sys- 
tem is 0.4e for noninteracting and 0.8e for interacting 
DW system exhibiting significant charge accumulation 
behavior. 

In Fig. 5, we show the interacting net charge current as 
a function of time for the case of magnetic field B = 0.5 T 
with different size of coupling window L w = (dotted 
black), 50 (dashed blue), and 100 nm (solid red). In the 
short-time response regime, shown in Fig. 5(a), the quan- 
tum interference dominates the time-dependent charge 
current feature with rapid oscillation. The oscillation 
amplitude and frequency of the time-dependent charge 
current remain similar for the cases with different win- 
dow size, this similarity is because the quantum inter- 
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FIG. 5. (Color online) The interacting net current Iqj versus 
time for L w = (dotted black), 50 (dashed blue), and 100 nm 
(solid red): (a) short-time response; (b) long-time response for 
B — 0.5 T. The other parameters are the same as Fig. 3. 



ference oscillation behavior is mainly due to the multi- 
ple scattering in the transport direction, and interference 
of subbands in the semi-infinite leads. In the transient 
switching-on regime t < 0.5 ps, the charge current for 
both the cases of short L w = 50 nm and long window 
L w = 100 nm are similar to the case without a window 
L w = 0.0 nm exhibiting the response time of the sys- 
tem from an isolated system to an open system. Later 
on, the charge current for the case of short window is 
suppressed by 0.4 nA, while the charge current is en- 
hanced for the case of long window by 0.8 nA. It should 
be noted that this quantitative feature is different when 
the e-e interaction effect is ignored, in which the charge 
current is almost the same for the cases without window 
L w — 0.0 nm and long window L w = 100 nm, however 
the charge current is suppressed by 1 nA for the case of 
short window L w = 50 nm (not shown). 

In the long-time response regime, shown in Fig. 5(b), 
the time-dependent charge current displays slow oscilla- 
tions and approaches to a steady current within 0.1 — 
0.2 nA. It is shown that the steady current is enhanced 
for the case of long window L w = 100 nm due to the 
Coulomb interaction. However, the Coulomb interaction 




FIG. 6. (Color online) The spatial distribution of the inter- 
acting many-electron charge density with different coupling 
window: L w — 0.0 nm (left), L w — 50 nm (middle), and 
L w = 100 nm (right) at time t = 14, 31.2, 100, and 200 ps. 
B — 0.5 T, and the other parameters are the same as Fig. 3. 



for the case of short widow is not significant on the time- 
dependent charge current in comparison with the pure fi- 
nite length DW system without window coupling. When 
the Coulomb interaction is ignored, the steady currents 
of the short and the long window are both suppressed 
(not shown). This demonstrates again the dynamics of 
the time-dependent charge current in the long-time re- 
sponse regime is significantly affected by the Coulomb 
interaction. 

To investigate how the window size affects the trans- 
port dynamics, in Fig. 6 we present the spatial distribu- 
tion of the many-electron charge at t = 14, 31.2, 100, 
and 200 ps, labeled by a-d in Fig. 5, respectively. It is 
clearly seen that, for both short and long window, the 
electrons perform inter-wire backward scattering in the 
short-time response regime (say, t = 14 and 31.2 ps), 
while the electrons are allowed to perform inter-wire for- 
ward scattering in the long-time response regime (say, 
t = 100 and 200 ps). This means that the former quan- 
tum interference dominant short-time response regime, 
the electrons favor the inter-wire backward scattering; 
while the latter Coulomb interaction dominant long-time 
response regime, the electrons favor inter- wire forward 
scattering. The many-electron charge density is mono- 
tonically increased in time. Furthermore, it is demon- 
strated that increasing window size can enhance not only 
the inter- wire scattering feature, but also the local charge 
accumulation at the coupling window. 
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IV. CONCLUDING REMARKS 

To conclude, we have performed a numerical calcula- 
tion of the time-dependent electric current and spatial 
charge distribution through a window-coupled parallel 
double quantum wire system based on GQME formalism 
including the electron-electron Coulomb interaction with 
the "exact diagonalization" method, and without resort- 
ing to the commonly used Markovian approximation. We 
have analyzed transient currents and their dependence 
on various parameters of the system with a certain ini- 
tial configuration and time-dependent switching-on cou- 
pling to the leads. For a given coupling window, we have 
demonstrated time-dependent transport properties of the 
noninteracting and the interacting DW systems. Apply- 
ing an appropriate magnetic field, we have found a short- 
time response regime dominated by quantum interfer- 
ence and inter-wire backward scattering. Moreover, the 
Coulomb correlation is significantly enhanced in the long- 



time response regime, 48 and the inter- wire forward scat- 
tering through the coupling window dominates the dy- 
namical transport properties. The conceived mesoscale 
window-coupled DW system could serve as an elementary 
quantum device for sensitive spectroscopy tools for elec- 
trons and quantum information processing by controlling 
the coupling window and the applied magnetic field. 
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